Correlated errors can lead to better performance of quantum codes 
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A formulation for evaluating the performance of quantum error correcting codes for a general error model 
is presented. In this formulation, the correlation between errors is quantified by a Hamiltonian description of 
the noise process. We classify correlated errors using the system-bath interaction: local versus nonlocal and 
two-body versus many-body interactions. In particular, we consider Calderbank-Shor-Steane codes and observe 
a better performance in the presence of correlated errors depending on the timing of the error recovery. We also 
find this timing to be an important factor in the design of a coding system for achieving higher fidelities. 
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I. INTRODUCTION 

The theory of fault-tolerant quantum computation (FTQC) 
has been developed as a realistic extension of the ideal theory 
of quantum computation. In an ideal prototype, all compu- 
tational components are functioning with no deficiency, while 
in the real world deviations from the ideal design is inevitable. 
This can be due to imperfections in the components or inter- 
actions imposed by the environment. In FTQC, schemes of 
error correction are combined with computational elements in 
a special design to prevent the spread of errors through a cir- 
cuit. Applying error correction steps in a multi-layer circuit is 
a well established procedure for FTQC |0 • 

Inspired by classical coding theory, a fault-tolerant form of 
quantum error correction codes (QECCs) can be achieved by 
concatenating single units of encoding. The performance of 
fault-tolerant codes has been studied for various noise mod- 
els. The classical model of a Markovian-independent noise 
has been thoroughly investigated in the literature 
In this picture, each individual component of the circuit is only 
coupled to a single designated bath, with no cross talk be- 
tween the baths. In addition, each bath is assumed to be large 
enough to constitute a local Markovian decoherence channel. 
Other possible defects coming from imperfections in compu- 
tational gates are also assumed to have local and instantaneous 
effects. All these give rise to a probabilistic description of the 
noise process. Upper bounding the error probability guaran- 
tees an arbitrary small likelihood of computational failure. 

However, an analysis of quantum computer proposals re- 
veals the inaccuracy of the above noise model. For exam- 
ple, ion qubits in an ion trap setup are collectively coupled 
to their vibrational modes |9[] . In a quantum dot design, dif- 
ferent qubits are coupled to the same lattice, thus interacting 
with a common phonon bath II 1 QTI . The exchange interaction is 
the main candidate for implementing two-qubit gates in solid- 
state proposals. Recently, it has been shown that many-body 
exchange interactions can be strong enough to act as a major 
source of noise [11]. Equally important, collective control of 
the qubits may also give rise to errors due to the inaccuracy of 
the control field [12]. These examples invalidate the assump- 
tion of error independence (locality in space). In addition, the 
non-Markovian nature of noise has been observed in various 



systems 111 311 . Therefore, the assumptions of exact locality in 
time and space should be relaxed to attain a more realistic 
model lfl4ll . First attempts to introduce correlations into noise 
models were in a classical manner: multilocation joint prob- 
abilities that are stronger than the independent model yj, |3[] . 
Recently a physical approach to the problem was taken by in- 
troducing a Hamiltonian description of the noise process 1 6, 7]. 
However, because of the discrete nature of error correction, 
the noise Hamiltonian was perturbatively treated to achieve a 
faulty paths description of the noisy computation. In another 
recent paper [8], the authors explicitly show the importance of 
a non-perturbative calculation of the FTQC threshold. 

In this paper, we study the performance of quantum codes 
in a single block of error correction towards a comprehensive 
analysis of the FTQC problem with a realistic model. Our aim 
is to compare spatial correlations in noise without any limiting 
assumptions on time. Focusing on well-defined Calderbank- 
Shor-Steane (CSS) codes [15], we first formulate a measure 
for the code performance and then apply it to different deco- 
herence scenarios. 

Our main finding is that a coding system can function bet- 
ter in the presence of error correlations. This is mainly due 
to the entangling power of a common bath which can help in 
the preservation of the code-state coherence designated for re- 
cov ery II 1611 . This is in contrast to results from previous studies 



II. CORRELATED ERRORS 

First we define the notion of correlation in errors. From 
the physical point of view, the primary probabilistic descrip- 
tion of the noise originates from two main assumptions on the 
system-bath interaction: each qubit is coupled to its own bath, 
and the resulting decoherence is described by a completely 
positive (CP) map. A CP-map modeling of the noise process 
allows a probabilistic interpretation of the errors which may 
be otherwise impossible 12111 . Based on this classical random- 
ness picture, correlations could be defined by some joint prob- 
ability distribution for the multi-qubit error operators. How- 
ever, a microscopic description is lacking in such a model, 
motivating a more physical formulation. We remark that the 
source of error is either a second system with a large Hilbert 
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space or an imperfect driving field of a quantum gate. 

In general there are three possible scenarios as depicted in 
Fig- El (0 errors with no correlation, which corresponds to 
model of separated baths, each coupled only to a single qubit; 
(ii) short-range correlations in which all the qubits are only 
coupled to a common bath; and (iii) The most general case 
with long-range correlations where qubits are interacting with 
each other while sharing a common environment. 

In the following, we locate the qubits of a code in the above 
scenarios and compare their performance as a function of de- 
coherence time Td- This is when we enter the recovery phase. 
The choice of CSS code is made to obtain an analytical solu- 
tion to the problem. A [n, k, d] CSS code is an n— qubit code, 
encoding k bits of information, capable of correcting errors 
acting on at most t = [^^-] different qubits. Consider an n- 
qubit quantum code in a single step of error correction. The 
n qubits are sent through a noise channel M expressed as a 
completely positive (CP) map: N{p s ) —J2 aEapsE^. We 
expand this map in the basis {£<*!,...,«„ = C3> ... <g> cr™^}, 



(a) 



1 Environment 
• Qubit 
Coupling 



* * * * 



(b) 



(c) 



Figure 1: A sketch of (a) an uncorrelated error model, (b) a short- 
range correlation, and (c) a long-range correlation. 



{p}=pi,---,p n , 

{g}=qi ,...,«„ 

where {do = I, <7i = a x , o-i = a y , 03 = <r z }. Follow- 
ing the matrix product structure of the basis, we introduce an 
wth-level submap M w (p s ) to be the sum of the terms in the 
above expansion (Q3 of which the array (pi + qi, ...,p n + q n ) 
has w non-zero elements. The matrix representation of the 
superoperator JV W = [e| u p | r g J has the following elements: 



— / e {p},{q} (Pi + Qi) nas w nonzero elements, 







otherwise. 



(2) 

The noise map AT(p s ) can be decomposed as a sum of sub- 
maps J\f w 

w w,{p},{q} 

(3) 



Theorem 1 Consider a CP noise map $> CP (p) = 
5^ F a pF^ and a code space Pc- A recovery CP map 
1Z correcting the map can also correct a Hermitian 

map = J2 a c a F a pFl. 

Proof. This is a straightforward application of theorem 1 in 
Ref. Q. ■ 

It is obvious that the correctable sum noise map JV C is a 
Hermitian linear map by definition. Therefore one can find a 
recovery map IZ(p) inverting J\f c , which belongs to the cor- 
rectibility domain of a [n, k, d] code,; i.e. there exists a real 
number px such that 



1Z o J\f c — pj^T. 



where X is the identity superoperator. 



(5) 



We can separate the correctable (by means of a [n, k, d] 
CSS code) and uncorrectable parts of the map by splitting 
the noise map into submaps Af c (ps) —J2 w <f^ w (Ps) m ^ 

Kc{ps) =Et< w ^(ps)- 

M(ps)=M c (ps)+M ic (ps)- (4) 

The correctibility of the linear map J\f c follows from the 
fact that a Hermitian linear map can be corrected by quan- 
tum codes j24ll . By a Hermitian map we mean a linear map 
$ : A4 U — ► M. v (where M. m is the space of m x to matri- 
ces) with a symmetric representation Q(p) = J2a c aFapFl, 
where the "map operators" {F a } are v x u matrices and the 
c Q s are real numbers. 



III. PERFORMANCE MEASURE 

We aim to study the performance of a quantum code for 
different decoherence processes; therefore, we need a proper 
measure to quantify it. As we discussed, the noise map Af is 
partially correctable by the recovery map 1Z. We introduce p^ 
in Eq. © to be the "performance measure" of a code. Before 
a rigorous derivation of this measure, let us apply it to a model 
of independent single-qubit noise channels 

M mdp {p s ) = {£ 1 ®...®£ n }{p s ). (6) 

where £. L (p) = (1 — p)p + pe(p), with p the probability of the 
error channel e. Suppose 1Z is the recovery channel correcting 
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errors in up to t different locations in the code; then we have 

K o J\P ndp (p s ) = no{Af^( Ps ) + NT c dp { Ps )} 

t 

= ^(l-p) (1 - c V^° U® ri ®e-e®^® rc }(/O s ) 



c=0 



n+...+r 

—n—c 



c=t+l 



(l-p) (1_c V^° E U 0ri ^e-e^^ rc }(/o s ) 

ri + ...+r c 

—n — c 

= ( E (") ( J - P) (1 ~V J + (error). 



By definition of the performance measure, p^i n d P is equal 
to Y^ c =o (™) (1 — p)^ 1_c ^ c - This value can be interpreted in 
the language of probability as the likelihood of having fewer 
than t individual errors M22I1 . In the above definition we dis- 
carded 7Z o Mic as the error. Actually this term vanishes in an 
averaging procedure introduced in Ref. 

10. 

Consider a code 

word ps — \^) (W| experiencing the noise Af(ps) and being 
corrected by a CSS code recovery map 



K(p) 



E 



-PS{ v }S{ A( }pS{ /1 }S{„}P. (7) 



jG{0,l}",Me{0,3} r 
\v\,\fj,\<t 



Our method is similar to the one introduced in a pertinent re- 
cent work IU7I1 for CSS codes, but we generalize it to an arbi- 
trary noise model. We choose fidelity 112211 as a proper measure 
of code performance: 



F(p s ,K(M(ps))) ={*\K(Ar(Ps))\*) 



(8) 



An analytical criterion is achieved by averaging the fidelity 
over all code words A CSS quantum code \C) is spanned 
by the code words | Q) where Q is a coset of C% and C^, (C^ 
is the orthogonal space to d, C2 C C\ C Z£). Applying 
the average over all pairs C2 C G\ suppresses all terms with 

e {ph{p} or e M#M which ^ of the order of 2 ~° {n) - We 
have verified this fact for a general noise model, justifying the 

bipartitioning of the noise map in Eq. (|4j. As a consequence, 

the fidelity average converges to a simple form of 



PX=F ave {\C){C\,1l{M{\C){C\)))= J2 



J {p},{p}- 



(9) 



n^i,{p} 



Such a code-independent measure equips us with a quanti- 
tative tool to compare different noise models. Suppose system 

and bath start from an initial product state pss(O) = ps(0) ® 
ps(0), with the bath density matrix pb{0) = J2i ^i\bi){bi\. 
For a system-bath propagator U Td (jd is decoherence time 
before recovery operation), the dynamical Kraus operators 

Now we can 



{Eij} can be written as 



^i{bj\U Td \h 



OAr— 
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Figure 2: Performance difference of a [7,1,3] CSS code for two sce- 
narios of local and non-local environments as a function of decoher- 
ence time r d . A negative sign p l ^f cal {gr d ) — p r ff nlocal {gr d ) demon- 
strates better performance for the nonlocal error.(iV = {7, 196}, 

pn = o.oi) 



find a closed form for p^ as a function of U Td and ps (0): 
- iE K\Tr^ {v} (b 3 \U Td \h))\ 2 



\v\<t 



PM(r d ) = ^ ^s(UlZ {v} )Tr s (Z {v} U Td )p B (0)}. 

\v\<t 

(10) 

In the following we apply this measure on different system- 
bath interaction configurations. 



IV. LOCAL vs NONLOCAL ENVIRONMENTS 

Here we compare system-bath configurations (a) and (b) in 
Fig. [TJ Each qubit is coupled to its own bath versus all qubits 
coupled to the same bath. We consider a spin star model which 
consists of a central spin-i particle representing a qubit, sur- 
rounded by N localized spin-i particles acting as a spin bath 
for the qubit B23I1 . In the local model the central spin a is 
coupled to the bath spins <7j, (i = 1, N) via a dephasing 
interaction of the form 



N 



Hhf = 2_j 9i a zZii 



(ID 



where a z and Z% are the z components of a and respec- 
tively. This is the effective Hamiltonian of dephasing in NMR 
due to spin-spin interactions 112511 . Also this hyperfine contact 
type of interaction is a source of electron spin decoherence in 
quantum dots due to the interaction with nuclei [26]. In the 
presence of a strong magnetic field along the z direction, the 
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dephasing interaction is the dominant term in the hyperfine 
interaction 027I1 . 

Two different configurations of local and non- 
local decoherence for an n-qubit code are mod- 
eled as Hi oca i = YT m =\TH=\9im.vf Z T and 

H non i OC ai = J2m=iJ2i=i 9im<rT z i- As an example 
we consider a [7,1,3] code surrounded by N = {7,196} 
bath spins [23]. The bath state, ps in Eq.dTob. is a thermal 
state of the spin bath at temperature T — lmK with an 
internal Hamiltonian Hb = 2j=i ^ z %- In addition, for 
simplicity we assume a symmetric coupling g im = g. The 
difference of the performance of local and nonlocal models 
P% cal {gTd) - P™ nlocal {9 T d) is plotted in Fig. We observe 
that the performance strongly depends on the decoherence 
time Td- Consequently, choosing the correct time for recovery 
can improve the functioning of the code in the presence of 
correlations. This is in contrast to a classical system in which 
a common source of error acts simultaneously on different 
parts, and therefore leads to a higher risk of computational 
failure, while in the quantum case, indirect coupling of the 
qubits through a shared common bath or direct coupling by 
means of many-body interactions imposes additional dynam- 
ics on the code coherence that changes the performance of the 
code. Furthermore, these results reveal that engineering the 
qubits of a code block in the same location (non-local bath) 
may have an advantage over spatially separating them (local 
bath), as suggested in Q7[| . 

Another interesting comparison can be made between lo- 
cal and nonlocal baths of the same size, with the correspond- 
ing Hamiltonians H local = J2m=i Ei=f 9im<*™ z T and 
Hnoniocai = YZn=i Si=i 9im<7? ' z i- The simulation results 
are plotted in Fig. |2] which demonstrates the time-dependent, 
alternative behavior of the code performance difference. 



V. TWO-BODY vs MANY-BODY INTERACTIONS 

We now introduce many-body interactions as a counterpart 
to multivariate joint probability in classical descriptions of 
correlations. In a recent study |7j], the fault-tolerant quantum 
computation in the presence of many body interactions has 
been investigated. Following their formulation, we include 
additional terms which act simultaneously on qubit pairs 
< ij > and on the bath. As an example consider the noise 
Hamiltonian 
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Figure 3: A [7,1,3] code performance in the presence of two- 
body and three-body interactions. The performance difference 

two — body { \ three — body { \ j . ■ , .< j?r * c *i 

PjV \ Td ) ~Pm \ Td > determines how the eirect oi three- 

body interaction varies by decoherence time Td (N — 7). 



Fig. [3] for short recovery times they behave similarly, while at 
longer times an alternative performance of two cases emerges. 

VI. FAULTY GATE 

A faulty gate is another main source of errors in compu- 
tation. In many proposals for quantum computers, a mag- 
netic field or laser pulse is an ingredient of the control field 
which can act locally [22] or globally (nonlocally) 11211 . Ob- 
viously an inaccuracy in manipulating the control field causes 
inaccuracy in the computation. Here we compare the fault 
tolerability of an n-qubit code driven by an imperfect mag- 
netic field acting locally (H m i = g Y17=i ^i z i) or glob- 
ally (H mg = gBo Y^i—i z i)- Because of the classical na- 
ture of the error, we assume a random process for the con- 
trol error, i.e. Bi — Bideai + W%, where the random 
variable Wi has a probability distribution pwi{w) = p(w). 
For a time r r the total rotation generated by this field is 
Umi{w,T r ) — exp[— ir r g Y17=i Bi{w)Zi] for local control, 
and U mg (w,T r ) = exp[~iT r gB (w)Y / " = iZ i \ for global 
control. The distance of two unitary operations U and V, 
acting on a Hilbert space 7i, can be measured by the m- 
erage fidelity Ave W€H (ip\V^U\ip) = ±\Tr(V^U)\ 2 , with 
d = dim(7i). To quantitatively compare a noisy gate with 
the ideal one, we use this distance averaged over the random 
variable Wi. For local control this value becomes 



Hi 



-body 



N 



9 2^2^ a i Zi 

J=l i=l 



j,k=l i=l 



in which g and g' determines the strength of two-body and 
three-body interactions. The ratio g' /g, in some specific con- 
figurations of spin-^ particles, can rise up to 16% ll ill . We 
compare two cases of two-body only interactions, with g' = 0, 
and three-body interactions with g' = O.lg. As shown in 



riocai( T rg) = ( / cos(r r guj)p(uj)dw) n , (13) 

J — oo 

while for the global one it is 



(cos ( T r guj ) ) n p ( 10 ) dhJ . 



Applying the Holder inequality 12811 . we find that 

riocal{ T rg) < Fglobal{ T rg)- 



(14) 



(15) 
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This indicates the higher fidelity of global control in compar- 
ison to the local one. 

VII. CONCLUSION 

In summary, we have studied quantum error correction in 
the presence of correlated errors. A microscopic descrip- 
tion of the noise process enables us to classify correlations 
based on their physical relevance in quantum computer de- 
signs. Roughly speaking, we find that it is not necessarily true 
that correlations in errors imply additional flaws in compu- 
tation, but on the contrary, these correlations may positively 
affect the performance. Furthermore, the results presented in 
this paper emphasize that, besides choosing a proper code, 
the pre-recovery decoherence time is a main factor to achiev- 
ing higher fidelity, since the code performance can drastically 
vary with the timing of the recovery operation. 
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